Superconducting-normal interface propagation speed in superconducting samples 
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In this paper a new approach to obtain the interface propagation speed in superconductors by 
means of a variational method is introduced. The results of the approach proposed coincide with 
the numerical simulations. The hyperbolic differential équations are introduced as an extension of 
the model in order to take into account delay effects in the front propagation due to the pinning. 
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I. INTRODUCTION 

The study of interface propagation is one of the most 
fundamental problems in nonequilibrium physics. The 
understanding of the magnetic field pénétration or ex- 
pulsion in Superconducting samples has been a major 
challenge. An important problem to be solved is the dé- 
termination of the speed at which the interface moves 
from a superconducting to a normal région. 

In RefQ], Di Bartolo and Dorsey have obtained the 
front speed by using heuristic methods such as Marginal 
stability hypothesis(MSH) and Réduction order. 

In gênerai, the nonlinear équations have been employed 
to model fronts propagation in différent areas such as 
population growth and chemical reactions. Our start 
point is the nonlinear diffusion equations(ND) of the form 
Ut = u xx + f(u) obtained from the Ginzburg-Landau 
expressions 2 (GL). The GL comprise a coupled équations 
for the density of superconducting électrons and the local 
magnetic field. 

Benguria and Depassier^£ have developed a varia- 
tional speed sélection method(BD) to compute the front 
speed in ND équations. In the BD method a trail func- 
tion g (x) is defined a priori and one may find accurate 
lower and upper bounds for the speed c. Only if the lower 
and the upper bounds coincide, then the value of c can 
be determined without any uncertainty. To eliminate the 
ambiguity in the speed détermination, Vincent and Fort 
m RefH have proposed a more accurate approach based 
on the BD method. The approach assumes some approx- 
imative considérations from where the function g (x) is 
determined. 

The purpose of this paper is to develop further the in- 
sights into the front propagation afforded by the work in 
RefQ. We aboard the détermination of the propagation 
speed from a variational point of view. We use an alter- 
native approach to the one developed by Vincent-Fort^. 

In order to describe the évolution of the System be- 
tween two homogeneous steady state, we assume a 
SC sample embedded in a stationary applied mag- 
netic field equal to the critical H = H c . The mag- 
netic field is rapidly removed, so the unstable normal- 
superconducting planar interface propagates toward the 
normal phase so as to expel any trapped magnetic flux, 



leaving the sample in Meissner state. We have considered 
that the interface remains planar during ail the process. 

The existence of a delay time in the interface propaga- 
tion Systems is an important aspect that can be modeled 
by hyperbolic diffusion equations(HD) which generalize 
the ND. The HD has been recently applied in biophysics 
to model the spread of humansl, bistable Systems^, for- 
est firesâ and in population dynamicsi^. With the goal to 
take into account the delay effect on the interface prop- 
agation speed 2 in superconductors, due to, for example, 
imperfections, vortex-vortex interactions, the présence of 
pinnin g 11 ! 12 , we have included the relaxation time r for 
the front, and indee introduce the hyperbolic differential 
équations. 



Traveling wave solutions. In this paper, we are inter- 
ested in the one-dimensional time-dependent Ginzburg- 
Landau équations, which in dimensionless units 2 - are: 
Ôtf = (l/« 2 ) d 2 J -q 2 f + f- f and âÔ t q = Ô 2 x q - f 2 q. 



Here, the quantity / is the magnitude of the supercon- 
ducting order parameter, q is the gauge-invariant vec- 
tor potential (such that h = d t q is the magnetic field), 
g is the dimensionless normal state conductivity (the 
ratio of the order parameter diffusion constant to the 
magnetic field diffusion constant) and k is a parameter 
which détermines the type of superconducting material; 
k < 1/ \/2 describes what are known as type-I supercon- 
ductors, while k > 1/V2 describes what are known as 
type-II superconductors. 

We are interested in finding traveling wave solutions 
for our model. We will search for steady traveling waves 
solutions for the GL équations of the form f(x, t) — s(x — 
et) and q(x, t) — n(x — et), where z — x — et with c > 0. 
Then the équations become 



—^s zz +cs z ~ns + s~s =0, 



acn. 



s 2 n = 0, 



(1) 



2 



II. VARIATIONAL ANALYSIS 

Vector potential q — 0. In this section, we assume 
= for the GL équations, 



dtf 



1 



dlf + f-f 



(2) 



Then, there exists a front / = s(x — ci) joining / = 1, 
the state corresponding to the whole superconducting 
phase to / = the state corresponding to the normal 
phase. Both states may be connected by a traveling front 
with speed c. The front satisfies the boundary conditions 
linis^-oo / = 1, linis^oo / = 0. Then Eq.((2]) can be writ- 
ten as, 

Szz+CK 2 s z +$ k (s) = 0, (3) 

where gfc = k 2 s(1 — s 2 ) and # = (I/k 2 )^. 

We definc p(s) = —ds/dz > and g such that 
h = —dg/ds > 0. Taking into account hp + {gfik/p) > 
2 \fgh$k and following the BD method- we arrive to 



c> {2/k) / (gh$)* ds/ 
Jo 



g ds. 



(4) 



Now, the asymptotic speed of the front for sufficiently 
localized initial conditions may be determined in the limit 
t — * oo. In the limit one has s —> \ for z — > — oo, and s is a 
slowly varying fonction of z. Thcrefore one has s zz <C s z , 
and from Eq.® we have that k 2 cs z + k 2 $(s) ~ and 
p ~ -s z ~ #/c. 

Assuming p = ^(s) /a > 0, where a is a positive con- 
stant to détermine, we can write in gênerai form the trial 
function as, 



g {s) = exp 



S" 1 ^) ds 



(5) 



Multiplying in both sides by the function h in the ex- 
pression ^Sug/p = hp, we have that, 



h$kg = h 2 p 2 



(G) 



By using Eq.©, the relation h^^g = h 2 $ 2 /a 2 is ob- 
tained. Then, the following relation is valid, 



2 \Jh 5 fc 5 = — h $. 



(7) 



Substituting Eq.([7]) in Eq.Q, the gênerai expression for 
the speed is given by, 

c~ maxa^o.j) (-^- J $(s)h(s) ds/ J g (s) ds^j .(8) 



Taking into account the form of $ and Eq.©, the trial 
function can be written as 



9(s) = [(s 2 1)/ 



,21« 2 /2 



(9) 
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FIG. 1: Illustration of the front speed obtained by différent 
methods in the q = case versus the GL parameter. 



and the function h(s), 

h(s) =a 2 S -( 1+Q2 )(l- S 2 )(" 2 - 2 )/ 2 . 
The intégrais in Eq. (J5J) are given by, 

f g( s )ds=-^T((l- a 2 )/2)r(l + a 2 /2) 
Jo V 71 " 

which is valid for a > 0, 

1 f(s)h(s)ds = ^[T((l-a 2 )/2) 



(10) 



(H) 



(12) 



2r((3-a 2 )/2)]r(l + c 2 /2), 

we arrive to the 

2T [1(3- 



which is valid for < a < 1. 

Replacing Eqs.fTT)) and (QU in Eq 
speed for the front, 



1 - 



k a 



(1 




(13) 



Notice that for a = 1, we obtain the maximum for 
Eq . (fT3]) , c — 2/k which is the resuit obtained by using 
the MSH method. 

In Fig.l the front speed versus the time delay is shown. 
The continuous line represents the results of the approach 
proposed in this paper following Eq. (fl~3|) and the numer- 
ical simulation by Eq.©. The results coincides. Also, 
the dashed line represents the bound from the varia- 
tional(BD) method^. 

Vector potential q = 1 — /. For a set of parameters^ 
k = \/\[2 and â — 1/2, we have that s(z) + n(z) = 1, 
then Eq.ÇO) takes the form s zz + (c/2) s z + #(s) = 0, 
where 3"(s) = s 2 (l — s) is the reaction term. Proceeding 
as in Eq.© we have that, 



g(s) = [(1 - s) /sf exp (a 2 /s) 
and the velocity is given by, 
'2 V2 



(14) 



max a£ i 



3(s)h(s) ds/ f' g(s) ds) .(15) 
Jo I 



3 



The interface speed is given by, 



c ~ max 



q£(0,1) 2 V2 



r(<* 2 ) 
r(i + a 2 ) ; ' 



(16) 



for a = 1/2 we obtain the maximum for Eq. (|16p . then 
c = \/2 which is the resuit obtained by using the MSH 
method. 



III. FRONT FLUX EXPULSION WITH DELAY 

It is wcll known the existence of pinning produces a de- 
lay timeii in the magnetic field pénétration o expulsion. 
This can be taken into account by resorting to hyperbolic 
differential équations seen in Section I, which generalize 
the parabolic équation. The aim of this section is to 
study of the interface speed problem in superconducting 
samples by means of the HD équations, which can be 
written as 



d 2 u du 



dt 2 dt 9: 



u , te \ , d f( u ) 



(17) 



In the absence of a delay time (r = 0) , this reduces to 
the classical équation u t = u xx + f(u). 

Vector potential q = 0. Taking into account the Eqs.Q 
and (|17|) we can write the following expression, 



K 2 T 



d 2 f 



*df Ô 2 f a 2 dS 

K "77T" — ~ — ~ ~t~ K tf + K (18) 



dt' 



dt 2 dt dx 2 

where # = s(l — s 2 ). 

It has been prove d 8 ' 9 ' 10 that Eq. (fï7f has traveling wave 
fronts with profile s(x — ct) and moving with speed c > 0. 
Then we can write Eq. (fT51) as follows, 

(l-ac 2 )s zz + c[ K 2 -ad\s)}s z +dk{s) = 0, (19) 

where z = x — et, a — k 2 t, 3fc = k 2 and with boundary 
conditions iim z _,ooS = 0, lim z ^.- 00 s — 1, and s z < in 
(0, 1); s z vanishes for z — > ±00. 

We define p(s) = —ds/dz > and g such that /i = 
—dg/ds > 0. Taking into account (1— a c 2 )hp+(g$k/p) > 
2Vl — ac 2 y/g and following the BD method we ar- 
rive to 



vr 



> 2k- 



a c 



Jo g{ k 2 - oï') rfs 



(20) 



In order to obtain the trial fonction g(s), we take in 
considération that in the lim s — > 1 we get 



c[K 2 -a3'(s)}p + S K (s)~0, 



(21) 



since s zz -C s z . Then, we write an expression for p in 
terms of 
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FIG. 2: Time-delayed interface propagation speed for the case 
of q — versus the time delay r. 



The expression for the speed is given by, 

2k 2 [h$/( K 2 -a3>)] ds 

Jo g( k2 - a ^0 ds 



1 — a c 2 



max ae ( ,i) 



(23) 



where the intégrais can be only solved by numerical meth- 
ods. Taking into account Eq. (|22|) and the expression 8 
(1 — ac 2 )hp — $ K g/p, we have obtained the relation for 
the trial fonction, 



g = exp 



(1 - ac 2 ) J $ K 
from where we have for our case, 



• ds 



(24) 



g = exp 
and for h 



2k 2 (1 -ac 2 ) 



(3as) 2 + lg 



( s 2 _ -Q(2a+K 2 ) 2 
„2(a-K 2 ) 2 



(25) 



h(s, a, a) 



a 2 (K 2 + a(3s 2 -l)) 2 
K 2 s(s 2 -l)(ac 2 -l) 



g(s,a,K) (26) 



In Fig.2 the front speed versus the time delay is shown. 
The continuous line represents the resuit of the approach 
proposed in this paper following Eq. (|23p which coincides 
with the numerical simulation done using Eq. (|18[) . Also, 
we have represented the lower and upper bounds from 
the variational(BD) method^. 

Vector potential q = 1 — /. Taking into account the 
Eqs. fTT]) and (fT5)) we can write the following expression, 



r d 2 f 1 df = 8 2 f 
2 dt 2 2 dt dx 2 



2*2 dt' 



(27) 



where 5 = s 2 (l — s). 

Then we can write Eq. ([27|) as follows, 



(1 



a c 2 ) s zz + c[k 2 ~ ad' (s)} s z + = 0, (28) 



P = 3 K /ot [k 2 - a$'(s)}. 



(22) where we have assumed §7: = (1/2)^ and a = r/2. 
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FIG. 3: Time-delayed interface propagation speed for q 
1 — / versus the time delay. 



The expression for the velocity is given by 



> 2V2-5! 



J L g(l-2aZ>) ds 



(29) 



In order to obtain an expression for the trial fonction 
g(s), we take in considération that in the lim s — > 1 we 
get 



:[{l/2) -a$'(s)]p + 5 K (s)~0, 



(30) 



since s zz -C s z . Then, we write an expression for p in 
terms of 5, 

p = $Ja[{l/2)-a$]. (31) 
The expression for the speed is given by, 

c 2 k 2 JoMftg/qi/^-aff')] ds 
- ~ max a£ ( d ; : (32) 



Taking into account Eq. pTj) and the expression 8 (1 — 
ac 2 )hp — ^ K g/p, we have obtained the relation for the 
trial function, 



g = exp 



((1/2) -a^ 



/\2 



(1 - ac 2 ) J 3 K 
from where we have for our case, 



■ ds 



(33) 



;(s, a, a) = exp 



4 s (1 - ac 2 ) 



(,9i +92) 



(34) 



where g\ = 1 + 24 a 2 s 2 (3s — 2), and g 2 — lg (s 

^sCla+l) 2 s s(16a-l) 



/i(s, a, a) 



a 2 [1 + 4as(3s- 2)] 
4s 2 (l-s)(l-ac 2 ) 



g(s,a,a) (35) 



The intégrais in Eq. ([32]) can be only solved by numerical 
mcthods. 

In Fig.3 the front speed versus the time delay is 
shown. The continuous line represents the results based 
on the approach proposed in this paper following Eq. |32|) 
which coincidcs with the numerical simulation done using 
Eq. (|27p . Also, we have represented the lower and upper 
bounds from the variational(BD) methocU^. 

Conclusion. We have computed for the Ginzburg- 
Landau équations in the form of parabolic and hyperbolic 
équations the superconducting-normal interface propaga- 
tion speed by a new approach. This approach is based 
in the method proposed by Vincent and Fort in Refd. 
We have obtained the expressions for the trial function 
g in each case developcd. The results of our methodol- 
ogy coincide with the numerical results for the examplcs 
analyzed. 
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